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We demonstrate the possibility for explicit construction in a discrete Hamiltonian model of an 
exact solution of the form exp( — |n|), i.e., a discrete peakon. These discrete analogs of the well- 
known, continuum peakons of the Camassa-Holm equation [Phys. Rev. Lett. 71, 1661 (1993)] are 
found in a model different from their continuum siblings and from earlier studies in the discrete 
setting [Phys. Rev. Lett. 83, 248 (1999)]. Namely, we observe discrete peakons in Klein-Gordon- 
type and nonlinear Schrödinger-type chains with long-range interactions. The interesting linear 
stability differences between these two chains are examined numerically and illustrated analytically. 
Additionally, inter-site centered peakons are also obtained in explicit form and their stability is 
studied. We also prove the global well-posedness for the discrete Klein-Gordon equation, show the 
instability of the peakon solution, and the possibility of a formation of a breathmg peakon. 
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I. INTRODUCTION 



In the last two decades, intrinsic localized modes (ILMs), also termed discrete breathers (DBs), have become a 
topic of intense theoretical and experimental investigation; see, e.g., [l| for a number of recent reviews on the topic. 
Per their inherent ability to bottleneck and potentially transport the energy in a coherent fashion, such exponentially 
localized in spacc and periodic in time entitics have comc to be of interest in a variety of contexts. These range from 
nonlinear optics and arrays of waveguides 2] to Bose-Einstein condensates (BECs) inside optical lattice potentials [Š| 
ö ' and from prototypical models of nonlinear springs |J] to Josephson junetions 5] and dynamical models of the DNA 
| double strand |fj. 

i The ubiquitous nature of nonlinear lattice systems (i.e., arrays of coupled nonlinear oscillators) has prompted the 
examination of the behavior of nonlinear waves in these systems, particularly for the waves that are well-known and 
understood in the continuum analogs of these equations i.e., in nonlinear partial differential equations. In this direction 
of aetivity, many of the coherent, nonlinear wave struetures that have previously been discovered in continuum settings 
such as regular solitons 0, compaetons JŠ[ , shock waves or gap solitons 0] have been recently examined in the 
discrete setting; see, e.g., the works of [J for regular solitary waves, f° r discrete compaetons, ^2] f° r discrete 
| shock waves or ^ŠJ for lattice gap solitons. 

On the other hand, a continuum nonlinear wave that was recently obtained in the context of shallow water wave 
equations, namely the peakon (a peaked soliton solution with a discontinuity in the first derivative at its peak) only 
very recently started to be considered in the discrete context |l4|. Peakons weregiven their name (in the context 
of the so-called Camassa-Holm equation, which is a dispersive, integrable model 0,0) due to their discontinuous 
derivative at peak amplitude. It is worth noting, however, that such "sharply pointed" struetures were also obtained 
much earlier in the context of nonlocal models for plasmas (see, e.g., 17]). Whilc, at first glance, it may not appear 
particularly natural to have discrete analogs of these solutions, as derivatives are not, strictly speaking, defined in 
the context of the spatial lattice, our aim in the present work is to examine the "discrete peakon" . This will be a 
lattice profile in the form (in fact, exactly) ~ exp(— a\n\) that in the continuum limit will asymptotically approach 
the continuum peakon solution. 

A number of twists accompany this discrete peakon solution introduced here. In the diseussion, we will illustrate 
that glimpses of solutions that can be eategorized under this new "species" may have already appeared in a somewhat 
generalized form in earlier works. We hope this may initiate a more general examination of the presently suggested 
lattice peakons. We believe that we have encapsulated the key mechanisms whose interplay can give rise to the 
existence of peakons (in both the discrete and the continuum setting) in a model of the Klein-Gordon (KG) variety 
and its nonlinear-Schrödinger (NLS) type analog. It is interesting to note that both the dispersive and the nonlinear 
terms that we have combined in the present setting have appeared previously in diverse contexts (that will be 
mentioned below) ; however, they were never combined to allow for the existence of peakons. Another intriguing bit 
concerns the stability properties of the peakons which are also examined in what follows. We find that the peakons, 
while unstable in the KG variant of our model, become stable in its NLS version. This is because the negative energy 
direction present the former model is prohibited by the charge conservation in the latter model. In the KG case, 
discrete breathing peakons are possible instead. Interestingly, such solutions were also identified earlier in elosed form 
in a elass of models very different than the ones examined herein [more specifically in a elass of homogeneous, nearest 
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neighbor, Klein-Gordon Hamiltonians] in pj . Finally, it is remarkable that the presently examined class of models 
and its stability features can be treated by methods available for the continuous nonlinear neid equations as illustrated 
below. 

Our presentation will be structured as follows. In Section [nj we will discuss the model and its motivation. In 
Section UIH the numerical results will be presented for site-centered, as well as for inter-site centered peakons. In 
Section IlVl we will examine the stability of such structures through analytical considerations based on constraincd 
energy minimization, as well as on functional analytic arguments. In Section the existence of discrete breathing 
peakons is discussed. Finally, in Section IVII we will summarize our findings and present our conclusions, as well as 
some open questions for future studies. 



II. MODEL AND MOTIVATION 



Examining the peakon from the "reverse engineering" (or the inverse problem) point of view, we would like to 
use the properties of a solution of the form u(x) ~ exp(— |x|), to construct a (continuum as well as a) discrete 
lattice with dynamics that supports such peaked solutions. Adopting this viewpoint, some of the key properties of 
u{x) = exp(— \x\) are that 

u'{x) = — u(x)sgn(x), (1) 
(l-d 2 x )u(x) = 25(x)u(x). (2) 

This second property (cf. also [lfij ) justifi.es why a strongly localized impurity (of the form of a <5-function) may be a 
relevant context in which such peaked solutions could arise. We will return to this point in Section fVTl 

Another, perhaps more interesting for our purposes, property is the result of convolution of such a peaked function 
with a Kac-Baker exponential interaction kernel J(\x — y\) = exp(— \x — y\) [2ÕI l2lj| . The convolution yields: 

/oo 
exp(-\x - y\)exp(-\y\)dy = (1 + \x\) exp (-\x\) . (3) 
-00 

This suggests immediately from a mathematical perspective a Klein-Gordon (KG), as well as a nonlinear Schrödinger 
(NLS) model with long-range interactions that would support exact peakon solutions of the form u{x) = Aexp(— a\x\). 
In particular the KG model would read: 



f°° ( 1 ( u 2 

u tt = aj exp(-a|x - y\)u(y)dy - I 1 - - ln I — 



(4) 



The peakons in this case would represent static solutions of the KG equation. The corresponding NLS model would 
be of the form 

f°° 1 /M 2 \ 

iu t = -a J exp(-a|x - y\)u{y)dy - - ln I — j u, (5) 

wherein the peakons correspond to standing wave solutions of the form u(x) — Aexp(it) exp(— a|a;|). 

Continuing along this lane of reverse construction (the motivation for each of the relevant terms will be given 
below), we now explain that an interesting feature of the above considerations is that they can also be carried through 
in the discrete setting. In particular, for u n — exp(— a\n\), we surn up the geometric series establishing that: 

E, , ,x / exp(2a) + 1 A 
exp(-a n - m \)u m = — — - + n \\u n . (6) 
me z Vexp(2a)-1 ) 

Consequently, for the discrete Kac-Baker interaction kernel 

J nm = exp(-|n - m\), (7) 

we can devise a discrete KG, as well as a discrete NLS model that have discrete analogs of peaked solutions. The 
discrete KG model is of the form: 
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exp(2a) + 1 1 j (u. 
exp (2a) - 1 ~ 2~a ^ [Ä 2 



(8) 
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with the exact discrete peakon solution 7r„ = Aexp(— a|n|), while the corresponding discrete NLS chain can be 
formulated as: 



I1i n — / u nm ^71 



Ih 



exp(2a) - 1 2a \ A 2 



U n , (9) 



where a discrete peakon given by the standing wave exp(it)ir n is the exact solution of the model. 

It is these latter equations [Eqs. (|SJ) and JÜJ] and their dependence on the parameter a that determines the interaction 
"range" that we plan on investigating in what follows. Notice that a can be considered as a natural spacing parameter. 

It is interesting to note as a side remark (to which we will return in later sections) that this model not only 
supports an exact "on-site" discrete peakon solution such as the one given above, but additionally supports exact 
"inter-site" peakon solutions of the form: u n — Bexp(— \n — 1/2|) for the Klein Gordon (in the DNLS case it is 
u n — B exp(ii) exp(— \n — 1/2|)). The value of the prefactor B is given by the relation 



\n(B/A) = a 



(exp(a) - 1) 



2(exp(a) + 1) 



(10) 



Such explicit solutions (especially inter-site ones) are rarely available in non-integrable discrete models. Inter-site 
solutions are typically less stable than their on-site siblings [l^- In the present setting, we will study in detail the 
behavior of such two-site peakons numerically as well as analytically. 

In motivating the model, aside from its intrinsic mathcmatical interest due to the existence of the peaked solutions 
(both in the discrete case and in the continuum limit), we should remark that both the dispersive and the nonlinear 
terms included here have appeared in a variety of settings before. The Kac-Baker type kernel [2ÕL l2l| , aside from its 
relevance in models of statistical physics, has been used quite extensively in recent nonlinear studies of lattice models 
emulating biopolymer dynamics including DNA; see, e.g., |22t l23l l24l I2ŠL l2o . |27 ] . Hence, this type of interaction is 
rather ubiquitous and can be controllably tuned (depending on the value of a), to be practically nearest neighbor 
(for large a) or much longer range (for a — > 0). Let us note also that in the past and in the framework of continuum 
equations similar forms of this kernel had been examined in the work of [2^ (but in a rather different dynamical model, 
namely one of the KdV type) wherein it was found that traveling waves acquired a peaked waveform. In contexts 
more closely related to the ones of the present work, let us also mention that "cusp solitons" (i.e., peakons) were also 
found in continuum models of the NLS type with Kac-Baker interactions in |29| (however, they were unstable) and 
in the case of a nonlocal Klein-Gordon field theory in [ŠÕj . 

The logarithmic nonlinearity in nonlinear model equations was originally introduced in the context of quantum 
field theory [3l| . It reappeared in |32j| . where it was proposed as an equation for generalized quantum mechanics. 
Later in |33| . it was suggested as a description for extended objects in nuclear matter, while more recently it was 
examined in the context of a scalar field model in inflationary cosmology j.llj . These studies have also triggered a 
more mathematically oriented interest in this nonlinearity and the properties of the solutions of the corresponding 
nonlinear wave equations [Š^- Perhaps, most closely to the purposes of interest to this study, this type of logarithmic 
field theory has appeared in saturable nonlinear optical media. The initial investigations of [36| in the latter context, 
in the framework of the "mighty morphing" spatial solitons, were later placed in a more physically realistic framework 
in connection with photorefractive materials in 37] . The work of [37J suggests that for the nonlinear waveguidc 
evolution, the logarithmic nonlinearity provides an accessible model that offers valuable insight, while maintaining 
the characteristic features of the underlying physical process. A note of caution should however be made in this 
connection in that the nonlinearity of Eqs. ©-10 should be viewed as a more reasonable physical model for larger 
amplitudes (where it can be considered as an approximation of a more physical nonlinear term such as ln(l + |u n | 2 )). 
For amplitudes tending to 0, the divergence of ln(|it„| 2 ) appears to be somewhat unphysical and leads to the absence 
of a small amplitude excitation (so-called "phonon" ) spectrum. 

The combination of the features of the dispersive interaction (its controllable range and wide applicability) and of 
the logarithmic nonlinearity (an accessible one representing adequately a number of physical processes) renders our 
model a possibly good playground to study, e.g., an array of coupled saturable nonlinear (logarithmic) waveguides. 
Both the potential relevance of our results in this context, as well as their inherent mathematical interest in establishing 
the discrete properties and behavior of the peaked solutions, lead us to examine Eqs. © and JjJJ in what follows. 
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III. NUMERICAL RESULTS 
A. General Setup 

The equations that we will examine can be re-written in a more general form: 

Ü n - JnmUm + F(u n ) = (11) 

mGZ 

for the KG lattice and 

iü n = - X! J nmU m + G{u n , u* n ) (12) 
for the DNLS chain; recall that J nm is given by Eq. 0. For Eqs. © and (JÜJ, the respective on-site terms are: 



F(u n ) 
G(u n ,u*) = 



exp(2a) + 1 _ J_ ^ fui 



exp(2a) - 1 2a \A 2 

2 1 / u n u* 

— In 



exp(2a) - 1 2a \ A 2 



(13) 
(14) 



Without loss of generality, we set A — 1. The exact solutions of interest for Eq. are of the familiar peakon form 
mentioned previously: 

7r„ = exp(-a|n|). 

We examine the linear stability of these solutions by using in Eq. i|ll|) 

u n = ir n + eexp(iu)t)v n , (15) 

where ir n is the original peakon and u> are the eigenfrequencies of linearization around the solution (v n are the 
corresponding eigenvectors). The resulting linear stability equation (obtained by using the ansatz of Eq. (|15fl to O(e) 
in Eq. (fTT|l ) reads: 



■u; 2 



V n = ^ J ™™ V ™ ~ F' (lT n )v n . (16) 



tit 



This is an eigenvalue problem for the matrix J nm — õ nm F' ' {~K n ). The discrete peakon is linear ly unstable if there are 
eigenfrequencies lo with the negative imaginary part. Since the matrix elements J nm are bounded and translation- 
invariant (only depend on n—m) while F'(Tr n ) exponentially decays as n — > oo, one can show that the eigenvalues of the 
truncated matrix, with \m\, \n\ < N, will tend to the eigenvalues of (jpp|l as N — > oo. The eigenvalue for the truncated 
matrix can easily be solved using numerical linear algebra packages; this gives the approximate eigenfrequencies lv 
and the corresponding eigenvectors v n . 

For the DNLS lattice, the stability can be performed in the "co-rotating" frame using the ansatz 

u n = exp(it) [ir n + e {a n exp(— iujt) + 6* exp(ia>*£))] . (17) 

Then, the resulting linear stability equations will be of the form: 

where J is the linear stability (Jacobian) matrix of the form 

duj du*. \ 
\ õu j du* J 

and T n = — Ylmez Jnm u m + G(u n , u n ) + u n (the Jacobian should be evaluated at the peakon profile, u n — 7r„). 
We now proceed to examine stability and dynamics properties of peakons in Klein-Gordon and DNLS systems. 
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FIG. 1: Left panel: Spatial profile of an exact discrete peakon (a = 1). 

Right panel: Spectral plane (Re(w),Im(w)) of the stability matrix for this solution in a Klein-Gordon chain. 



B. Klein— Gordon 1-site peakons 

The profile of a peakon is shown in Fig. Q together with the spectrum of eigenfrequencies of the equation linearized 
at the peakon. For the KG chain, we find that the solutions are unstable (for all values of a). This is because of a 
negative energy direction that leads to an imaginary pair of eigenfrequencies. The corresponding eigenvector has the 
same shape as the peakon itself. 

This fact can be observed in more detail in Fig. [21 where the dependence of the imaginary part of the unstable 
eigenvalue of the stability matrix is shown as a funetion of a. This figure also gives the dependence of the energy 
of the peakon on a, which can be analytically calculated: E(a) = A 2 coth(a)/(2a). From these figures, it can be 
dedueed that the solution becomes less unstable with (inereasing) a, or, in other words, when the width of the peakon 
decreases. This can be equivalently interpreted as a weaker instability as the solution approaches its anti-continuum 
limit, single-site peakon sibling. This is a rather natural feature of spatial disereteness which typically serves to 
stabilize coherent struetures that are unstable in the continuum limit (e.g., due to collapse) or even ones that do not 
exist in that limit 0. 

In order to examine the dynamical evolution of the instability of the KG lattice peakon, we used direct numerical 
simulations, performed with a 5th order Calvo's symplectic integrator |38| with a time step At = 0.001, which 
preserves the energy up to a faetor 10" 15 . We introduce a perturbation = eir n (to the exact peakon solution n n ) 
with e = 0.1 to excite the unstable eigendirection. The exponential growth of the peakon instability is shown in the 
left panel of Fig. [3] If the perturbation were = — e-K n with e = 0.1 again, the peakon does not grow. Instead, it 
evolves to a breathing state (see e.g. right panel of Fig. |3J), which will be analyzed in SectionlVl 



C. DNLS 1-site peakons 

As explained above, the (spatial dependence of the) profile of a DNLS peakon is the same as that of its Klein-Gordon 
analog (see Fig. However, in the DNLS setting, the peakon is stable for all values of a. This fact has its origin 
in the gauge invariance of the solutions of DNLS-type equations. In particular, an interesting feature of the discrete 
peakons is that the U(l) symmetry of the DNLS chain prohibits the single negative energy direction that was present 
in the KG lattice. Essentially, the unstable direction of the KG lattice is transversal to the same-charge hypersurface 
in the DNLS case. As a result, perturbations along this potentially unstable direction are banned by the presence of 
the extra symmetry. 

Figure 01 shows the spectral plane for a typical case together with the dependence on a of the charge (also referred 
to as power in optics) of the peakon. The charge is defined as Q(u) — J2 n \ u n\ 2 /2 and it can be observed that its 
value decreases with a and tends to Q = 1/2 (as should be expected as a — ► oo). This dependence can be analytically 



6 




'0 2 4 6 8 100 2 4 6 8 10 



a a 

FIG. 2: Left panel: Dependence of the maximum imaginary eigenfrequency (i.e., maximum real eigenvalue) of the stability 
matrix on a. 

Right panel: Dependence of the energy of the peakon on a. 




FIG. 3: Left panel: Instability evolution of a discrete peakon at different times (a = 1). Exponential growth can be observed, 
while the energy of the system is preserved up to a 10~ 15 precision. The different snapshots of the solution correspond (from 
inner to outer) to times 0, 1.8, 2.7 and 3.6. 

Right panel: Time evolution of a perturbed peakon that develops into a breathing state. The displacement of the central sites 
of the peakon are shown as a function of time. 



calculated: Q(a) = Ä 2 coth(a)/2. 

The stability of the solution can be verified in the time evolution numerical experiment shown in Fig. \5\ which has 
been performed through a 4th order Runge-Kutta integrator with time step At = 0.01. The phase space plot at the 
central site shows that a randomly perturbed solution remains orbitally close to the exact discrete peakon solution. 
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FIG. 4: Left panel: Spectral plane of the stability matrix for the peakon of Fig.^in the case of the DNLS chain. 
Right panel: Dependence of the peakon charge (power) Q — 7r^/2 as a function of a. 




Re(u(t)) Re(u(t)) 

FIG. 5: Left panel: Phase space diagram of the central site of a perturbed (full line) and an unperturbed (dashed line) DNLS 
peakon. 

Right panel: A blow-up of the left panel that illustrates the orbital stability of the peakon solution (since the perturbed solution 
remains in its vicinity). 



D. 2-site Peakons 



We now examine the behavior of two-site (i.e., inter-site) peakons, both in the KG, as well as in the DNLS chain. 
The energy and charge of 2-site peakons can be analytically calculated as E(a) — B 2 /(Aa sinh(a)) and Q(a) = 
B 2 /(2sinh(o)). 

Contrary to their single site counterparts, two-site solutions are unstable (also in the DNLS model). In this case, 
two negative energy directions and hence two imaginary eigenfrequencies could be identified in the spectral plane 
of the linearization eigenfrequencies in the case of the KG lattice, while one such eigenfrequency was present in the 
DNLS setting (see Fig. |BJ. The eigenmode corresponding to the KG case is antisymmetric. 

We also simulated the dynamical developmcnt of these instabilities, observing that KG two-site peakons are com- 
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FIG. 6: Spectral plane of a the two-site peakon (with a = 1). Left panel: KG peakon. Right panel: DNLS peakon. 
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FIG. 7: Time evolution of the two-site peakon in the DNLS chain (with a = 1) induced by perturbing the central particle. 
Left panel: The full line represents |uo| 2 , the dashed line represents |iti| 2 . 
Right panel: Charge density in space-time. 



pletely destroyed (as their one-site counterparts are also not stable). For the two-site DNLS peakons excited with the 
perturbation e<5„.o with e ~ 10~ 4 , the solution oscillates between the one-site and two-site peakons. These results are 
shown on Fig. [7| 



IV. ANALYTICAL RESULTS 



The above results motivate us to examine the stability of the Klein-Gordon and DNLS peakons from an analytical 
perspective and, in particular, using energetic considerations. We now proceed to study the stability of the uniform 
steady state ("vacuum") and of one-site peakons in each of these settings. The conclusions about stability or instability 
do not depend on the values of A and o, so we set 



A = l, o=l, 
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so that the peakon profile is given by 

7r„ = exp(-|n|). 

A. Klein-Gordon 

1. Hamiltonian formulation. 

We can rewrite the Klein-Gordon equation (|11J) as 

ü n + d Un T(u)+d Un W(u) = (18) 

where 

2»=-~ E e-l"- m l U „ Mm , (19) 

^) = E((t + \)<-\<*<)- ( 2 °) 

Above, A is a positive constant taken to be 

A = (.,.) = J] e - 2 H = J±j, (21) 



where (u, u) — ^2 n ^z u n u n (hi the Klein-Gordon case, we assume that the components u n are real-valued) . 

Remark 1 Later we will show that the Klein-Gordon equation llfy) is globally well-posed in Z 2 (Z) x Ž 2 (Z); if both u(0) 
and ü(0) belong to l 2 (Z), then there is a global solution u(t) with ||u(t)||ja < oo, ||ü(t)||/2 < oo for < t < oo. See 
Theorem^ At the same time, the norms ||it(t)||;2, ||ü(t)||ja could grow unboundedly large with time. 

We can rewrite (|18l) as 

ü + T'(u) + W'(u) = 0, (22) 

where T'(u), W'(u) may be interpreted as variational derivatives with respect to u of the functionals T and W. 
The value of the energy functional 

e kg (u,ü) = j2y +t{u) + w{u) (23) 

is conserved along the trajectories of (f22J - 

2. Stability of vacuum. 

First, let us notice that the zero solution is stable with respect to F-perturbations of the initial data. We bound 
T(u) by 



\T(u)\ < \ 



sup E e"!"-™! 



(«,«)< ^±y(«,u). (24) 



mGZ 

This inequality is due to the Schur test applied to the matrix J nm = e _ ' n_m ', which yields 

i i 

|| Ju\\i2 < i SUp \J nm \ ) ( SUp V" | J nm \ | ||u||ja = : r||u||p. (25) 
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In the expression (|20() for W(u), the first term is (^ + j){u, u). When the amplitude of u is small, the main contribution 
in l|20l) comes frorn the last term and is of positive sign; thus, for = sup m \u rn \ < e, 

Ekg(u) > (-t^tt + ^ + ]) («, «) + w ) ln 3. ( 26 ) 



2(e-l) 

which is positive for e sufficiently small. Thus, the zero solution minimizes the value of the energy functional among 
perturbations with bounded amplitude. Since 

1/2 

HIj.. = supK| < [ > >„| 2 J =|Hp, (27) 



sup|u„i < y2 \u n \ 2 



the zero solution also minimizes the energy among all perturbations with sufficiently small l -norm. 



3. Instability of one-site peakons. 



Now we address the stability of the peakon 7r n = e '"L Following Derrick |4l|, let us consider the family of vectors 
7r( A ) with the components 7Tb = A 1 / 2 ^, where n n is the peakon profile. We have: 



E kg (tt w ) = T(tt (a) ) + W(tt (x) ) =aA-6AlnA, 



with a = T(tt) + W(ir) and b = \ Enezjl > 0- 
Since 7r is a stationary solution to l|22|) . so that 



T'(tt) + T^'(tt) = 0, 



we have: 



d_ 

dX 



T(7r( A >) + W{n^) = (a- 6(lnA + 1)) 



r (A) 



a-& = 0, 



(28) 



(29) 



(30) 



dA 2 



T(7r( A )) + V^(7r( A )) = -- 
A 



= -6 < 0, 



(31) 



thus the stationary solution n does not minimize the energy. 

It is also easy to check directly that perturbations in the direction of the vector ^ | A _ 1 7r< - A ' 1 = \^ are linearly 
unstable. For this, let us consider the solution of the form 7r„ + r n (t); the linearized equation on r is 

r„ + ( T » + W"(7r)) nm r m = 0, 

where (T"(tt) + W"(ir)) nm = ^J^J M' Thc matrix r "( 7r ) + VF"(7r) has the eigenvalue -1, with the corresponding 
eigenvector being 7r itself: 

(T"(n) + W"(ir)) nm ir rn = -n n . 

The perturbation in this direction will grow exponentially, in accordance with our numerical results (cf. Figs. ITl3*|l . 
According to 42|, the linear instability in this model gives rise to the dynamic, or nonlinear, instability: 

Theorem 1 The stationary peakon solution u(t) = ir to ill'L'1) is unstable with respect to l 2 -perturbations of the initial 
data. 
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B. DNLS 



1. Hamiltonian formulation. 

We now turn to the DNLS equation l|12|l . which can be rewritten as: 

iÜ„ = 2d K E(u,u*), (32) 

where E(u, u*) — T(u, u*) + V(u, u*), with 

T(u,u*) = ~ £ e- |n - m| < Um , (33) 

(n,m)eZ 2 

v(u,u*) = E (f f - 1) I"™! 2 - il M »l 2ln W a ) ■ ( 34 ) 

nez V V / / 

The value of the energy functional 

E(u, u*) = T{u, u*) + V(u, u*) 
is conserved along the trajectories of (|32[1 . The value of the charge functional 

is also conserved (due to the t/(l)-invariance of the system). 



2. Stability of vacuum. 



Stability of the vacuum solution u = with respect to / 2 -perturbations of the initial data is proved in the same way 
as in the case of Klein-Gordon equation: The solution u = delivers the smallest (zero) value to the energy functional 
among all vectors ip of sufficiently small / 2 -norm. 



3. Stability of one-site DNLS peakons. 



In the case of DNLS, the perturbation u t— > u + eu would lead to smaller values of the energy functional as in the 
case of Klein-Gordon equation. However, now this perturbation is prohibited by the charge conservation. Therefore, 
we expect that the peakons are local minimizers of the energy functional under the charge constraint and hence are 
orbitally stable. 

We will prove the following theorem: 

Theorem 2 The standing wave solution e lt ir to (|12fl is orbitally stable with respect to l 2 -perturbations of the initial 
data. 

Let us remind the definition of the orbital stability (see 

Definition 1 The n-orbit {e IS 7r; s £ R} is stable if for all e > there exists õ > with the following property. If 
i*o G 1 2 (7j, C) is such that \\uo — 7r 1 1 ^2 < S and u(t) is a solution of Eq. \12^ in some interval [0, ii), then u(t) can be 
continued to a solution in < t < oo and 

sup inf \\u(t) - e ts 7r||;2 < e. 
o<t<oo s e R 

Otherwise the ir-orbit is called unstable. 
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Remark 2 Eq. ijJüj) is locally well-posed in 1 2 {7L, C) since /or u £ l 2 the right-hand side of $12\) (where we set A = 1, 
a = 1) is also in l : 



1 



ln {u n u* n ) 





< 







e+1 



e - 1 e 2 - 1 



+ MMH 



(35) 



We used inequalities $25}) and \2 r $ . Hence, for any uq £ Z 2 (Z,C), there exists r > so that there is a unique solution 
u{t) defined for < t < r that satisfies u(0) = uq. Due to the conservation of the charge Q{u) ~ ||u|| 2 2 /2 along the 
flow of Eq. 112]) . we conclude that u(t) is defined globally: u £ C' 1 ([0, oo), ^ 2 (Z, C)). This settles the question about 
the global well-posedness for $12\) in Z 2 (Z, C) which is a necessary condition for the orbital stability. 

According to we wm know the orbital stability of the peakon if we can prove that H = E"{-k) + Q"(ir) defines 
a quadratic form that is positive-definite on vectors that are tangent to the hypersurface of the same charge and are 
orthogonal to the orbit spanncd by 7r. 

Remark 3 Once one knows that H defines a positive-definite quadratic form on vectors tangent to same charge 
hypersurface and orthogonal to the orbit of tt, one proves that there exist 8 > and C > such that for any u with 
Q(u) = Q{tt) and \\u — 7r||p < S one has E(u) — E(tt) > C inf s ||u — e ls 7r||;2 (Theorem 3.4 in \lA l) and then the orbital 
stability follows from Theorem 3.5 in tl A l. 

We will use the real-valued formulation. For u(t) = v(t) + iw(t) £ ^ 2 (Z, C), with v = {v n }, w = {w n } real-valued, 
we write 



u(t) 



v(t) 
w(t) 



u(t) G l 2 (Z,R 2 ). 



The equation on u is 



ü = JE'{u) = 



o r 






-1 




_ V W E _ 



and the stationary equation on tv 



Let ip = £ + irj, ip = 



is given by 

E'(tv) + Q'(tt) = 0. 
The vectors tangent to the same charge hypersurface satisfy 



(36) 



while the vectors orthogonal to 
by e lt n) satisfy 



- Jtt (this vector corresponds to in, a tangent direction to the orbit spanncd 
(-Jtt,i/>) =-X) 7r » 7 ?»=0. (37) 



As we mentioned in Remark [21 Theorem |2 will be proved if we can show that the quadratic form defined by the 
Hamiltonian operator 



is positive-definite on vectors i\) 



H = E"{n) + Q"(n) 



, where both £ and r\ are orthogonal to tt. 



Let us find the explicit expression for H . For the second Frcchet derivative of T, we compute: 

(T"(7T)</>,V>>=- E ^ ml (Um+VnVm)- 

(n,m)£Z 2 



(38) 
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For the "potential" term, we compute 

(V"(n)^, tf) = £ (( A - + €) + \n\(C + Vi) ~ C) ■ (39) 



nt:; 



For the charge functional, we have 



(Q»^> + (40) 

Hence, the operator H = E"(-k) + Q"(tt) = T"(tt) + V"(tt) + Q"(tt) is given by 

<Ht/>, ^) = - E ( „, m)e z 2 e- |n - TO| Knfm + Mm) (41) 

+ E„ ez ((A-i + N)(^ + ^) + ^)- 



We have: 
where 



' H+ 




"C" 


#- 







(42) 



#+ m - -e-l"- m l + (A - 1 + |n|) £„ m (43) 

^ n - m =-e-l"- m l + (A+|n|)<5 nm . (44) 

Let us consider the symmetric matrix 

A nm = - e -l"-" l l + \n\6 nm + 7r„7r m . (45) 

Numerical computations show that this matrix does not have negative eigenvalues and that the dimension of its null 
space is dim N(A) = 3. One can check that N(Ä) is spanned by the vectors v^ + \ v^', and t>(°): 

=ö(n+l/2)e-l n l, 

'(-n+l/2)e -|n| , 



„(-)=„(+) 



= <W 

where 0(n) = 1 for n > and for n < 0. In particular, 

Tr = « (+) + - « (0) e iV(A). 

Since A is symmetric, the eigenvectors that correspond to other (positive) eigenvalues of A are orthogonal to tt. 
The matrices H+ m , H~ m can be expressed as 

H tm = A nm ~ 7T„7T m + (A - l)5 nm , (46) 
^nm — A nm TT n TT 7n ~t~ A$ nm — -\- S nm . (47) 

We have: H + tt = Air — (tt, tt)tt + (A — 1)tt = —tt, since (tt, tt) — A, and H ~tt = 0. Since H + and H~ are symmetric, 
their other eigenvectors are orthogonal to tt. But thcn they also have to be the eigenvectors of A. The corresponding 
eigenvalues are those of A shifted to the right by (tt, tt) — 1 > (for H + ) and by (tt, tt) (for H~), and hence are strictly 
positive. 

Remark 4 The important feature of this model that makes the analysis simple is that H + and H~ can be diagonalized 
simultaneously. 

Thus, assuming that ip = £,+ir/ £ i 2 (Z, C) satisfies (|36|l and H37[) (both £ and r/ are orthogonal to tt) and is different 
from zero, we conclude that (Hxjj, ip) is positive-defmite: 

(Hlp,1p) = ^ ( H nm^m + H~ m r) n r) m ) > 

(n,m)6Z 2 

and hence the one-site peakon solutions are dynamically stable. 
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4- Two-site DNLS peakons 



One can also prove the following theorem: 



Theorem 3 The 2-site standing wave solution to <|12|) is orbitally unstable with respect to l -perturbations of the 
initial data. 

The generalization of the matrices H + and H~ given above in Eq. (|43fl and 1441) is given by: 

H+ m = - e -l"-™l + (A - 1 - InK)) S nm (48) 
H- m = - e -l»-™l + (A + InK)) <W (49) 

In the case of the inter-site solution, known explicitly as 

Un = J^^e«e-l"- l /=l, (50) 
2 e + 1 

the eigenvalues of H + , H~ can be explicitly computed. As expected, u n is itself an eigenvector with eigenvalue — 1 
(see above), while the second eigenvalue is A2 = —0.149. However, since the eigenvalues of H~ are the ones of H + 
shifted by 1, this results in n(H + ) — n(H~) — 2, where n(H ± ) denotes the number of negative eigenvalues of the 
matrices . From the theory of |43L EL l45| (see also the recent work of ^B>E3l)j this implies that there is one real 
eigenvalue in the spectrum of linearization of the DNLS 2-site pcakon, in agreement with our numerical observations 
of Scction lTTTl This completes the proof of Theorem [3] 

V. BREATHING PEAKONS IN THE DISCRETE KLEIN-GORDON EQUATION 

A. Existence of Breathing Peakons 

We are interested in solutions to the Klein-Gordon equation <|18|) that have the form 

u„(t) = g(t)f n , (51) 

where g{t) is a scalar-valued function. 

Trying the peakon profile, f = n, so that u n (t) — g(t)e~\ n \, we get the following equation on g(t): 

9 = \g^9 2 : 9 = 9(t)- (52) 



We can integrate Eq. lO, getting 



^ + ^(l-lng 2 ) = £, (53) 



where £ > is the "energy of breathing". The function V(g) = -^-(1 — lng 2 ) (see Fig.|SJ) represents the potential in 
which g lives. 

B. Global well-posedness for the discrete Klein-Gordon equation in l 2 x l 2 

Before analyzing the stability of the breathing peakons, we need to consider the initial value problem for the 
equation of the Klein-Gordon lattice (fT5|l . 

Theorem 4 Eq. is globally well-posed in l 2 x l 2 . That is, for any initial data (uq^vq) € l 2 {1i) x Z 2 (Z), there is 
a unique solution u(t) that satisfies it(0) = uq, ü(0) — vq; this solution is defined for all times t > 0, and moreover 
||u(t)||;2 + ||ü(t)||ja remains finite for all t > 0: 

iieC 2 ([0,oo),l 2 (Z)). 
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FIG. 8: Potential V(g). The (unstable) stationary solutions g = ±1 correspond to the (unstable) peakons, u(t) = g(t)n = ±7T, 
while the (stable) stationary solution g = corresponds to the vacuum (which we know to be stable) . There are also oscillating 
solutions that correspond to breathing oscillations with the energy £ , < £ < 1/4. A breathing "bound state" at £ = 0.2 is 
also shown by dash-dotted line. 



Let us first discuss the local well-posedness of (|18|) in l 2 . We claim that for any (uq, vo) £ l 2 {^) X / 2 (Z), there exists 
r > such that there is a unique solution uit) defined for < t < r with (u,ü) = (ito,i>o) and u £ C 1 ([0, r), ^ 2 (Z)). 
This claim immediately follows from the observation that if u G / 2 (Z), then ü given by l|18fl is also from l 2 (Z). This 
is verified as in 

Now we turn to the global well-posedness. Let us prove the boundedness of L(t) = \\u(t)\\ 2 2 . Let us derive the 
equation for d 2 L(t)/dt 2 . We have: 



—L{t) = 2{ü,ü) + 2(u,ü). 

The total energy \{u 1 ü) + T(u) + W(u) is conserved along the flow generated by (|18|) allowing to express 

(ü, ü) = 2E - 2T(u) - 2W(u), 

where Eq — (vq,vo)/2 + T(uq) + W(uq) and (uq, vq) is the initial data that corresponds to the solution u(t). 
Thus, we obtain: 

^jL(t) = 4Eo-4:T{u)-W{u)-2(u,d a T)-2(u,d u W) 
= AE Q - 8T(u) + ( 4A + !) u l + 2«n ■ 



(54) 



Using the straightforward bounds on the terms in the right-hand side, namely 

e + 1 

that follows from 124(1 . and also 

where in the second inequality we used relation l|27|l . we conclude that, for some C = C(Eq) > 0, 

r/ 2 



dt 2 



L < C(l + L + 2LlnL). 



(55) 



(56) 



(57) 
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For L > 1, the right-hand side is monotonically increasing (and exceeds its range for < L < 1). Thcrefore, 
< L(t) < Z(t), where Z(t) is a function that satisfies 



and the initial data 



Z"(t) = C{1 + Z + 2Z\nZ) (58) 
Z(0) = max(l, 1/(0)) = max(l, (ito, uo)), 

= max (1, 2(uo, v )) . 



dL 

tf(0)=ma* 1,_ 



We can rewrite as 



Z" + C*ö z (-Z - Z 2 lnZ) = 0; (59) 
multiplying by Z' and integrating in t, we get 

CZ - CZ 2 \nZ = E, 



(Z') 2 



where E — ^ z ^ — h C(— Z(0) — Z 2 (0) ln Z(0)) is a constant of integration. Expressing Z' and separating variables, 
we get 

rZ(t) dz 



Z(0) 



VE + CZ + CZ 2 lnZ' 



Since the integral f^? n x , ^ z =g diverges at the upper limit (as ln 2 Z), we conclude that Z can not become 
infinite in finite time. 

The finiteness of Z 2 -norm of ii follows from relation 154fl . bounds (|55() and (|56() . and the finiteness of L = \\u\\ 2 2 that 
we already proved. 

This finishes the proof of Theorem^] 

C. Linearized stability of Breathing Peakons 

We will now analyze the linear stability of the breathing peakon, showing the absence of the exponential instability 
of small perturbations of the initial data. We rewrite Eq. (|18|) as the first order system, 

v = -d u (T(u) + W(u)), [W) 

and consider the perturbation of the solution (uo{t),vo(t)) — {g{t)ir 1 g{t)n) that corresponds to the peakon: 

u(t) = u (t + j(t)) + Su(t), v(t) = v {t + 7(t)) + Sv(t). 

The function j(t) adjusts the location of the breather (g(i)Tv, g(t)ir) so that it is closer (in a certain sense) to the 
perturbed solution (u(t), v(t)). We consider the linearization of system (JSÜJ. For this, we first compute 

T"(gn) + W"(gn) = T» + W"(n) - ^ = H+ - ^f, (61) 
where H + was introduced in (|46() . and then we can write 

g'(t + 7(*))7(*)t + 9täu(t) = öv(t), 



r+ :n„S, (62) 

We split 



g"(t + i{t))i(t)* + d t öv(t) = -{H+ - ^-)5u(t). 
öu(t) = a(t)n + <p(t), Sv(t) = b(t)ir + ip(t), 
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with (f>, ip £ 1 2 (1j) both orthogonal to tv. Projecting system 162(1 onto tv, gives the following system: 



g'(t + i(t))j(t)+ä{t) = b(t), 
g"(t + 1 (t))j(t)Tv + b(t) = -(-l 



(63) 



where we used the fact that tv is an eigenvector, H + tv = — 
the following system: 

/ 0(t)=^(t), 



Projection of i|62[l onto the direction normal to tv gives 



Üt) = -(H+ - 



(64) 



The analysis of system (I63f) is straightforward. We would arrive at this system if we pursued the stability analysis of 
Eq. (|52|l . At the same time, we could analyze that equation topologically. It corresponds to an unharmonic oscillator; 
its phase portrait in the (g,g) plane contains a set of closed trajectories circling around the origin (they correspond 
to the initial data (g,g) such that \g\ < 1 and the value of £ in (|53|l is smaller than 1/4). Each of these closed 
trajectories is stable with respect to small perturbations of the initial data: There is a neighboring closed trajectory 
passing through the point that corresponds to the perturbed initial data. 
Let us now analyze system (|64|l . We can rewrite it as 4> = JH<J>, where 



$ = 



1> 



1 

-1 



H 



H- 



_ ]ng_ 

2 

1 







with 4>, ip orthogonal to tv. The matrix H + (see l|46[) and thereafter) has eigenvalues 

<7(#+) = {-l>A-l,...}, 



where the only negative eigenvalue —1 corresponds to tv and the next eigenvalue, A— 1, is positive. The term — ^y- m 
(|61|l further shifts the spectrum upwards as long as \g\ < 1. Therefore, H is positive-definite on the space tv 1 - x ^ 2 (Z). 
Thus, \/H is well-defmed; the matrix JH is similar to \/H J\/H and hence has purely imaginary spectrum. 
This shows that the breathing peakon solutions g{t)Tv to (f 1 H|) are spectrally stable. 

Remark 5 This linearized approach to the stability does not prove the dynamic orbital stability of the breathing 
peakons. The nonlinear terms may transfer the energy between "breathing" oscillations in tv -direction and the per- 
turbations in the space tv , pumping the energy from system ffify) into |6'^| ). It is therefore possible that the energy 
of the breathing peakon, after a small perturbation, would wind up transferred, partially or completely, into directions 
orthogonal to tv (and then maybe back). We do not have a satisfactory description of this process, even though we 
believe that techniques such as the Hamiltonian dispersive normal forms of fA3 I could be relevant in addressing it. 
While outside the scope of the present study, this may be an interesting question for future investigations. 



D. Numerical Results 



An orbit of frequency lj v for g(t) can be determined by solving Eq. (|52|) . This can be done through a variety 
of methods such as, e.g., a shooting method in real space or using quadratures. Here we have chosen a Chebyshev 
quadrature method in order to integrate the equation. 

Figure [ü| shows the dependence of g(0) and £ of the peakon frequency ui p . It can clearly be observed (see the left 
end of the graphs) that as uj p — > 0, the energy approaches 1/4 and the amplitude is 1, hence the solution is very close 
to the unstable critical point of £ = 1/4 and g = 1. We can also observe that the peakon frequency can have any 
value as there do not exist resonances with the continuous spectrum (actually, such small amplitude, extended wave 
excitations do not exist since V"(0) = oo). 

Figure lTül shows the time evolution of a breathing peakon. It is worth pointing out that the main difference between 
DNLS peakons and KG breathing peakons is that, in the first case, only the first Fourier coefficient is different from 
zero, whereas in the second case, there are more non-zero coefhcients. 

We have confirmed in the numerical simulations that initial conditions corrcsponding to a perturbed discrete 
breathing peakon stays orbitally close to the exact solution. Figure ^] shows the difference between the evolution of 
the central partide of the peakon in a perturbed and an unperturbed case. The perturbation used is £ n = sõ n fi with 
e = 0.01. The perturbed peakon appears to be orbitally stable. 

Contrary to the static case, breathing two-site peakons are not destroyed by perturbations. Instead, the energy 
density oscillates as shown in Fig. ^| 
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FIG. 9: Dependence of g(0) (left) and £ with respect to the peakon frequency. 
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FIG. 10: Time evolution of the breathing peakon (the displacement of each of the first few sites is shown as a function of time). 



VI. DISCUSSION 

In this paper, we have engineered a mathematical model that has discrete pcakons as exact solutions. The continuum 
version of the model was also given and its ability to support the continuum analog of the solutions was highlighted. 
Both the dispersive and the nonlinear part of the relevant interactions were connected to earlier works. Furthermore, 
it was advocated that for suitable choice of the interaction range, this can be a relevant model in nonlinear optics 
with the characteristic features of photorefractive materials |37j| , while being a lattice dynamical model of interest in 
its own right. 

We have identified the discrete peakon solutions analogous to their (discontinuous in the first derivative) continuum 
limit, i.e., 7r„ ~ exp(— |n|), as well as their inter-site siblings. However, a natural question of interest would be whether 
there is a more general way of defining such solutions in the discrete setting. This is particularly relevant as solutions 
similar to the ones obtained here have appeared in other contexts. Such exam ple s consist of, e.g., the waveform of 
Fig. 1 in |40| (arising from the presence of an impurity) or that of Fig. 7 in |27j (arising because of the interplay 
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FIG. 11: Left panel: Phase space diagram of the central site of a perturbed (full line) and an unperturbed (dashed line) 

breathing peakon with to p = 1. 

Right panel: Blow-up of the left panel. 

The cross indicates the initial point of the simulation. 




FIG. 12: Time evolution of the two-site breathing peakon chain (with a = 1 and ui p — 1) perturbed with the antisymmetric 
mode. 

Left panel: The full line represents eo, the dashed line represents ei, where e n is the energy density of the n-th site. 
Right panel: Energy density in the space-time evolution. 



of nonlinearity with long range interactions, as is the case in this paper). Perhaps, an alternative criterion possibly 
involving the sign of the second difference close to the center of the wave could be used for a more general definition 
of the discrete peakon. This would be an interesting topic for future studies. 

We have also investigated the stability of such discrete waves and have found discrete peakons to be particularly 
interesting from this aspect as well. We have numerically observed that in the Klein-Gordon lattice setting such 
solutions are always unstable; however the negative energy direction responsible for this instability is eliminated due 
to an additional symmetry (the phase invariance that leads to the ž 2 -norm conservation) in the case of the DNLS 
chain. We showed (using a Derrick-type argument) that the peakon is not a local minimizer of the energy in the Klein- 
Gordon case and moreover indicated the direction of the perturbation that leads to a linear instability. We also showed 
that in the J7(l)-invariant DNLS equation the peakon is a local minimizer of the energy under the charge constraint 
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and hence is orbitally stable. Contrary to their one-site counterparts, two-site peakons proved to be unstable, as was 
explicitly demonstrated via a well-known functional analytic criterion rclcvant to DNLS type equations. 

Finally, using a separation of variables approach, we were able to show the existence of the exact periodic (breathing) 
peakon solutions in the Klein-Gordon lattice. The breathing peakons were shown to correspond to the subcritical 
initial conditions, while supercritical initial conditions lcad to the amplitude of the solution tending to infinity (in 
infinite time), with the unstable static peakon being the separatrix between the two types of behavior. We also 
systematically tackled the initial value problem for the Klein-Gordon lattice and showed that the finite-time blow-up 
of the F-norm of the solution is not possible, so that the system is globally well-posed. We then proved the absence 
of the linear instability of the breathing peakons. Yet, the question of the long-time behavior of perturbed breathing 
peakons remains open. 

It may be intcrcsting to try to extend this class of models to higher dimensional scttings and observe how their 
dynamical behavior is affected by the dimensionality of the underlying lattice. 
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